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Abstract 

A simple state sum model for fermions on a 1-manifold is con- 
structed. The model is independent of the triangulation and gives 
exactly the same partition function as the Dirac functional integral 
with zeta-function regularisation. Some implications for more realis- 
tic physical models are discussed. 

1 Introduction 

Discrete models are used in physics as a practical way of denning a functional 
integral. In lattice gauge theories one obtains a model that approximates the 
continuum theory better as the number of discrete lattice cells is increased. 
In the continuum limit, the gauge theory should of course regain the full 
rotational and translational symmetry. 

Topological quantum field theories can be denned using discrete models 
called state sum models [1] . In this case the partition function is independent 
of the number of discrete cells (the simplexes of a triangulation). Similar 
models are used in quantum gravity, but the requirement that the partition 
function is independent of the triangulation is generally relaxed; however 
the independence must at least hold in a suitable limit in which the model 
approximates general relativity. 
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Realistic quantum gravity models should of course include matter fields 
and so it is important to investigate models with matter fields that are in- 
dependent of the triangulation, either exactly, or in a limit. This is our 
motivation for studying state sum models with fermionic fields. 

In this paper a very simple one-dimensional state sum model is presented 
in which a fermionic field is coupled to a background gauge field. A sim- 
ple formula for the partition function of this model on a triangulated circle 
(i.e., a polygon) is presented in Section [2j It is demonstrated that the par- 
tition function is independent of the triangulation and depends only on the 
holonomy of the gauge field. 

A heuristic argument is presented in Subsection 13.11 that shows that the 
state sum model has an action that is a discretisation of the continuum 
Dirac action for a massless fermion field coupled to the gauge connection. 
Then, in Subsection 13.21 and the two appendices, the partition function of 
this functional integral is calculated precisely to show that it is in fact exactly 
the same as the partition function of the state sum model. 

The results are extended in Subsection 13.31 to models with a mass term. 
A state sum model with a discretised mass term is presented. However in 
this case the model does depend on the number of edges N in the triangu- 
lation of the circle. It is shown that as iV — > oo it converges to a limit, so 
that the triangulation- independence is recovered in this limit. The limiting 
value is closely related to the Dirac functional integral with a mass term, 
differing only by a certain phase factor. Finally, the paper ends in Section H] 
with a discussion of the results, including a comparison with operator cut-off 
regularisations of the Dirac functional integral. 

The results presented here complement previous work constructing quan- 
tum gravity state sum models with fermion fields in dimension three [21 E] 
and four (U [5] . These works construct discrete analogues of the continuum 
Dirac functional integral according to the heuristic continuum limit, as con- 
sidered here in Subsection 13. 1[ but do not have a direct comparison with 
the partition function of the continuum functional integral. The results pre- 
sented here give the first precise comparison of a discrete fermionic model 
with the partition function of the continuum functional integral. 

It is an interesting question as to whether our results can be generalised 
to a higher-dimensional model. Some related properties of the Dirac operator 
with a gauge field on a graph have been studied previously (6]; this suggests 
there may also be an extension of our state sum model to graphs. 
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2 The state sum model 



2.1 Fermionic variables 

In a functional integral, the fermion fields anticommute and so are elements 
of a Grassmann algebra. Therefore the state sum model is based on a finite- 
dimensional Grassmann algebra over C. Such an algebra is determined by a 
number of generators ai, a 2 , ■ ■ ■ , (fy satisfying the relations 

OjOj + OjOj = 0. (2.1) 

A function of the generators f(ai, . . . ,ai) is a polynomial, which termi- 
nates at the highest monomial a;a;_i . . . a\ due to the anticommutation. The 
Berezin integral 

dcti da 2 . . . da; / (2.2) 

is defined to be the coefficient of a«a;_i . . . a\ in the expansion of /. No 
independent meaning is attached to the differentials in this formula, and they 
do not appear outside an integral. However the order of them in the integral is 
important; transposing two neighbouring differentials in the notation changes 
the sign of the integral. 

It is possible to extend the above definitions to integration over a subset 
of coordinates, and perform the integral iteratively. So if / = a k ak-i . . . a±b, 
where 6 is a polynomial in the remaining variables a,k+i, ■ ■ ■ ,ai, then the 
integral is 

J dai da 2 . . . da k f = b, (2.3) 

with terms multiplying lower degree monomials in a± . . . a k integrating to 
zero. An example of iteration is the formula 

J dai da 2 f = J dai ^ da 2 f^j . (2.4) 

The integration enjoys translational invariance: if c is any polynomial in 
the generators excluding a±, then 

J /(ai + c) dai = J /(ai) da ± . (2.5) 



3 



In the applications of interest here the generators of the Grassmann alge- 
bra occur in pairs Oj, 6j and are also grouped into a number of n- dimensional 
vectors, such as 

ip = (ai,a 2 , ■ ■ .a n ), ip = (b x , b 2 , . . . , b n ). (2.6) 
In this case the integral is defined with the notation 

J dtp dip — J ddi d6i da 2 db 2 . . . da n db n . 

Let M be an n x n matrix with entries in C. The basic gaussian integral 

is 

J dip d^ = det M, (2.7) 

which is proved by expanding the exponential. 

The result (j2.7p can be extended to the case of fermionic source terms. 
Take c, d to be n-component vectors with Grassmann- valued entries that 
are polynomial of odd degree in the remaining generators (i.e., excluding 
components of both ip and ip), and M now an invertible matrix. We then 
have 

J dip dip e W+^+?d = det M e-civf" 1 ^ ( 2 8 ) 

which is proved by first completing the square with the translations 

V^^-cM" 1 , (2.9a) 
iP^iP- M' x d, (2.9b) 

and then using (12.71) . 

2.2 Definition of the state sum model 

The construction of the fermionic state sum model is based on fermionic 
variables at a discrete set of points. For this, a number of n-dimensional 
vectors ipi, ip i are used (with i labelling the vectors, not the components). 
The Grassmann algebra is the one generated by all of the components of all 
of the vectors. 

The basic construction is the partition function for the interval [0, 1] cR 
decorated with a fixed nxn matrix Q, as depicted in Figured! The partition 



4 



*« 



Figure 1: The partition function of an edge. The arrow indicates the orien- 
tation of the edge. 

function then is 

= e -^oWi_ (2.10) 

This has fermionic variables ip , associated to the boundary point 0, and "01 
associated to 1. 

Gluing two such partition functions together is carried out using the fol- 
lowing proposition, which states that one can multiply matrices by the use 
of Berezin integration. 

Proposition. 

r dij 1 d^ 1 e-^ Ql ^ 1 e i '^ 1 e-^ Q2i ' 2 = e -^ QlQ ^ 2 . (2.11) 
Proof. We have 

d^id^ e-^^e^e-^ 2 ^ 2 = J d^d^e^ 1 "^ 1 ^ 1 -^ 2 ^, 

(2.12) 

where I is the n x n identity matrix. The result (12.111) follows using (12. 8p 
with 

c = -%Q u d=-Q^ 2 , M = I. (2.13) 

Our interpretation of formula (12. lip is that there is a bilinear form on 
the fermionic states, 



(/, g) = J d^id^ /(^i) e^Wi), (2.14) 
and using this bilinear form to glue the partition functions results in 

Z [0*1] Z [0?1] = Z [0,l]' ( 2 - 15 ) 

with Q = QiQ 2 . 
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Clearly this can be iterated for the multiplication of any finite number of 
matrices, yielding the definition of the state sum model. Explicitly, 

J dipt . . . dV>;v-i d^ JV _ 1 e-^^e^e-^ 2 ^ 2 . . . e~^ Q ^ N 

= e -^ QlpN = Zf Q1] (2.16) 

with Q = QiQ 2 ■ ■ • Qn- 

The left-most expression in (I2.16P is to be interpreted as the definition 
of the fermionic state sum model on an oriented interval that is triangulated 
using N 1-simplexes. The equation itself is a statement that the partition 
function is independent of the triangulation, though of course provided the 
matrices Q and the orientation match. Geometrically, if the Qi are all invert- 
ible, they may be viewed as the parallel transport operators for a connection 
1-form on the interval. Then the partition function can be viewed as a func- 
tion of the connection that is invariant under gauge transformations in the 
interior of the interval. 

The state sum model on the left-hand side of (12.161) can be modified to 
include observables, that is, non-trivial functions of the intermediate variables 
ipi, In this sense the model is richer than the evaluation of the partition 
function on the right-hand side. The physical content of the model will be 
investigated further in Section [3] by comparing it to a functional integral. 

The partition function for the circle can be computed by gluing together 
the endpoints of the interval, as depicted in Figure HJ 




Figure 2: The state sum model on a triangulated circle. 

Mathematically, this is done by including an extra factor e^ N ^ N in the 
integrand, identifying ip = ip^, tp = ip N , and integrating over the newly 
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introduced variables, 



Z% = J d^ N d^ N e ^^ N e-^» QlpN = det (J - Q) , (2.17) 

where the last equality follows from ( 12.70 and the observation that VjvV'jv 
and ip N Qip N commute. 

An immediate consequence of (12.170 is that iBi vanishes if Q has an 
eigenvalue equal to 1, as is for example the case for the gauge group SO(2n + 
1). For SO(2n), Z^ is real- valued, as the eigenvalues occur in complex 
conjugate pairs. There are however gauge groups for which Zgj is complex- 
valued. An example of particular interest is U(l), for which we have Q = e~ ie 
with 9 G M, and so 

U(l) : Z% = l-e- ie . (2.18) 

This example shows that the partition function depends on the orientation 
of the circle for some gauge groups. Given a circle with a connection, if the 
holonomy around one direction is the matrix Q, the holonomy around the 
opposite direction is Q~ l . For SO(n), det (I — Q) — det (/ — Q^ 1 ) for all Q, 
and so the partition function is independent of orientation. For U(n), however 
det (/ — Q) is the complex conjugate of det (/ — Q^ 1 ), and the U(l) example 
(12.180 shows that the partition function is sensitive to the orientation for 
generic Q. 

2.3 Gauge transformations 

Gauge transformations act on the partition function of [0, 1] by a linear trans- 
formation acting on each set of fermionic variables. Thus if Uq, U\ are invert- 
ible matrices then the transformation is ip^ = U\^)\ and ip Q = (U )~ 1T iJ; , us- 
ing the inverse transpose matrix. The matrix Q transforms as Q' = UqQU^ ~ . 
It is clear then that the partition function for the interval (12.101) is invariant 
under gauge transformations, 

Similarly, the partition function of the circle (12.170 is gauge-invariant. 

Nothing so far has suggested that a variable ^ is related to tp^ However 
if, as the notation suggests, one regards these as complex conjugate vari- 
ables, then the corresponding transformations should be complex conjugate 
matrices, (Ui)~ 1T = Ui, which means that these matrices are unitary. 
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3 Comparison with the functional integral 



The partition function for the circle with a connection is compared to the cor- 
responding functional integral in this section. First it is shown, in a heuristic 
fashion, that increasing the number of 1-simplexes can be interpreted as the 
state sum integrand converging to the exponential of the Dirac action for a 
fermion field. Then a precise evaluation of the Dirac functional integral is 
compared to the state sum model. 



3.1 Continuum limit of the action 

This subsection demonstrates that the state sum model is a discrete version 
of the Dirac functional integral. This is done by looking at the action for the 
state sum model and taking the limit N — > oo, assuming that the discrete 
values of the field sample values of a differentiable field tp(t). The result of 
this heuristic limiting process is the continuum Dirac action. 

The state sum model for the partition function on a circle triangulated 
with N edges is 



= / II d ^ d ^ e^»*+e**..e*»-*»-^»-M»e+>>+». (3.1) 




The integrand can be rewritten as 

j — ^ e i>Ni>N e —tp N Qii>i^j ^ e i , ii>i e -i , iQ'2i>2 j ^ e i ! N-ii , N-i e -i> N -iQNi>N j 

= e ^UM^^), (3.2) 
with ip N = ipQ. 

With the choice At = l/N, where I is a positive constant, this looks 
precisely like a discretisation of the Dirac action on a circle of circumference /. 
Let t be the coordinate on the circle with range to I. Then the j-ih vertex 
corresponds to the coordinate t = jl/N, and At is the distance between 
neighbouring vertices. 

The limit N — > oo is the limit At — > 0. In this limit the parallel transport 
of the connection is Qj = 1 — iA a (t)X a At + O(At) 2 , where iX a are the 
generators of the relevant Lie algebra and A a (t) are the components of the 
connection 1-form. 

Assuming the field values are differentiable, we have 

^(ftH^M^), ,3,, 
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where the gauge covariant Dirac operator Ip is given by 



U) = ij t +A a {i)X a . (3.4) 

The single gamma matrix is equal to the complex number i. There is no spin 
connection contribution to ( 13. 4 j) because the tangent space is one- dimensional 
and the Lie algebra so(l) is trivial. 

In the limit At — > the sum converges to an integral, 

lim At^f(jAt)= / fdt, (3.5) 

giving the continuum limit of the state-sum integrand as 

hm J = jti^mm. (3. 6 ) 



3.2 The functional integral 

In this subsection the above heuristic considerations are confirmed by an 
explicit and precise evaluation of the functional integral in the XJ(n) case. 
The continuum partition function is 

Fgi = J V^V^e'ti^®**®, (3.7) 

where J]) is given by ( 13. 4 p and the fermions obey periodic boundary condi- 
tions. These boundary conditions can be regarded as the choice of the 'Lie 
group' spin structure for the circle [H]. The functional integral in ( 13. 7p is 
defined via zeta-function regularisation [7J [S], as reviewed in Appendix IA1 
The result is 

= det(0) = e ^ (0) e^ 2{0) , (3.8) 

where the functions rj^ and Cm 2 are defined in Appendix IA1 
It is shown in Appendix iBl that for U(l) 

F2i=z£, (3.9) 

where Q is the holonomy of the connection A. The result (13. 9p generalises 
immediately to U(n) by diagonalising the connection with a gauge trans- 
formation, whereupon the functional integral is the product of a number of 
U(l) functional integrals. The Dirac functional integral is invariant under 



9 



these gauge transformations because it depends only on the eigenvalues of 
the Dirac operator, which are gauge invariant. 

The result ( 13. 9 j) is surprising because the eigenvalues of the Dirac op- 
erator are unbounded and so the naive determinant of the Dirac operator, 
the product of its eigenvalues, diverges. Somehow the discrete model both 
approximates the eigenvalues of the continuum operator yet also avoids the 
divergence, and miraculously imitates the zeta-function regularisation. 

Some insight into the result ( 13. 9 p can be gained by comparing the eigen- 
values for the continuum Dirac operator with the eigenvalues for its discrete 
version. 

The discrete version of the Dirac operator is a matrix M acting on the 
vectors ip = ®f =1 ipj- It is determined by 



JV-l 



ipiMt/; = i>j Oj - Qj+iipj+i) 



3=0 



and can be written in block form as 



/ 1 



iM 



-Q2 
1 



\ 



-Qz 



1 - 



V -Qi 



Qn 

1 



(3.10) 



(3.11) 



where each Qi is in U(n). Since we have seen that the partition function 
associated with M is exactly the same as the continuum partition function, 
it must be that the matrix M is in some sense approximating the differential 
operator Ij). We now make this more explicit. 

For simplicity, we specialise to U(l). The eigenvalues of iM are fik = 
1 — a>k, where a& are the N roots of 



Q 



a 



N 



(3.12) 



and the corresponding eigenvectors are 



/ Qr 1 i \ 

ockQ 2 Qi 
<%Qz x Q2 X QT x 

\ a»- l Q- N l . . . Qr 1 J 



(3.13) 
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Taking Q = e~ l0 with 6 e [0, 2tt), as in (IZTSjl . we have 



' 8 + 2k-K 




fi k = l-e-^-^), (3.14) 

where the distinct eigenvalues are obtained by selecting a suitable set of 
distinct values of k. For example, this can be done explicitly by taking 
k = [(1 — N) /2], . . . , [(iV — l)/2], where [x] stands for the largest integer that 
is less than or equal to x. 

To compare with the continuum Dirac operator, we note that the eigen- 
values (I3.14p of %M have the large N expansion 



(3.15) 



while in Appendix [B] it is shown that the eigenvalues of ilp are 

il]): fi k = ii j, keZ. (3.16) 

The expressions (I3.15P and (13.161) coincide modulo 0(N~ 2 ) if the circle has 
length / = N and k is held fixed. The matrix iM hence approximates the 
operator ilp in the sense that the eigenvalues of small modulus coincide in 
the limit of a large circle. 

Note from (13. 15p that the eigenvalues of M are complex but the imaginary 
part is subdominant as N — > oo with fixed k. This raises the question as to 
whether it is fruitful to think of M as a cut-off version of Jp. We shall return 
to this question in Section HI 



3.3 Mass term 

In this subsection we show that inclusion of a mass term breaks triangulation 
independence, but gives a well-defined partition function as iV — > oo. We 
start now with a massive continuum fermionic action 

S= f dtijj(t)(]f>-m)ip(t), (3.17) 
Jo 

and we follow the usual discretisation procedure using (13. 3p and (13. 5p . The 
state sum is given by 

'JV-l \ 

Y[ dijjd^j I e Ef=o 1 ^((i-imA*)^-o i+1 v i+ i) ) (3 lg ) 
J=o J 
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where At = l/N, ip N = ip , and the mass term has been discretised according 
to 



o 



l __ N-l_ 

mi/jip dt -> mAt ^ ( 3 - 19 ) 

3=0 



We can evaluate Igi in a similar way as before to obtain 

I s i = det ^(1 - imAt) N - J| Q t j , (3.20) 

which is clearly not triangulation independent. The triangulation dependence 
appears to have crept in with the introduction of a fixed mass scale m. A 
way to regain triangulation independence might be to relegate m from a 
fixed parameter to a function of N and a new fixed positive parameter w! by 
1 — 2m At = e~ iml / N ; however, this would mean allowing m to be complex- 
valued. 

With m fixed, the limit N — > oo does however yield a well-defined answer, 
lim det |(1- imAtf - TT Q t | = det( e - imZ - Q), (3.21) 



i=l 



where Q is the holonomy around the circle. This can be compared to the 
evaluation of the partition function with a mass term in the continuum U(n) 
case. The calculation is done identically but using A' = A — m. The answer 
comes out as 

Wg[ = det(l — e 27Tim Q), (3.22) 

which differs from (I3.2ip by a phase factor. 

The disparity is due to the different way in which the gauge field and mass 
term are treated. In the continuum calculation they are bundled together 
into one term. However, in the discrete calculation, the gauge field is realised 
as parallel transporters while the mass term couples fields at the same point. 



4 Discussion 

In this paper, a fermionic state sum model on one- dimensional manifolds 
with a connection is defined. The partition function of the state sum model 
on the circle is the same as the functional integral of a complex fermion field 
with the Dirac action and the usual coupling to the connection gauge field. 
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A curious feature of the partition function for the circle is that it does 
not depend on the length of the circle, whereas the eigenvalues of both the 
continuum Dirac operator J}) and the discrete Dirac operator M clearly do. 

Some more insight into this can be gained by examining an operator cut- 
off regularisation of the functional integral. Such a regularisation depends 
on a cut-off scale c. It replaces iD with f(iD), where the function f(z) 
of a complex variable is the identity function f(z) = z for \z\ <C c, but 
f(z) = 1 for \z\ 3> c. Thus it effectively removes the eigenvalues of D with 
magnitude above the cut-off c. An example of such a cut-off regularisation 
is the Schwinger proper time regularisation. 

As c — > oo, the regularised determinant diverges. A calculation for the 
Schwinger proper time regularisation shows that the leading asymptotic term 
for log det \lp\ is (Z/7r)clogc. This can be confirmed in a simple way for 
special cases, e.g. A = 1/2, using a sharp cut-off and Stirling's formula for 
the asymptotic expansion of the factorial. 

The cut-off regularisation only agrees with the zeta function method once 
the leading asymptotic terms are removed [7]. It is worth noting that these 
leading divergent terms are proportional to I, thus explaining the role played 
by the length of the circle in evaluating this determinant. Thus to get a 
partition function that converges as c — > oo, one can multiply the cut-off 
regularised determinant with a 'cosmological term' e A ^ 1 , choosing a suitable 
function A(c) to renormalise the coefficient of I as c — > oo. This is the same as 
adding a term A(c)l to the exponent in (13. 7p when using an operator cut-off 
regularisation. The functional integral formula for this regularisation is then 



Similar divergences do not occur with cut-off regularisations of the eta invari- 
ant in the phase term [9] , so there are no additional parameters to renormalise 
for the phase of the partition function. 

Our conclusion from this discussion is that the state sum model is a more 
subtle regularisation of the determinant of the Dirac operator than a mere 
operator cut-off. 

Another aspect of this is the absence of the phenomenon of fermion dou- 
bling which commonly occurs with lattice discretisations of the Dirac action 
|10j . One of the standard assumptions of discussions of fermion doubling is 
that the action of the lattice theory is real (i.e. equal to its conjugate with a 
suitable definition of complex conjugation for fermionic variables). However 
the action used in the state sum model (I2.16P is in fact not real, and the 
resulting complex nature of the eigenvalues of M appears to be a crucial 
aspect of the model. 
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Finally, it is worth making a brief mention of possible generalisations. 
One can include an integration over the U(n) gauge field in the state sum 
model using the normalised Haar measure on the group 



The resulting partition function T[ 0j i] is independent of any background struc- 
ture; fT2TT5D yields 



so that T[o,i] is a projection operator. The corresponding partition function 
for a circle is thus an integer, and can be interpreted in terms of counting 
states in the space of polynomials in the fermion variables at a point that 
are invariant under the group action. 

Some aspects of the model presented here have been studied for the Dirac 
operator on a graph; the formula f)3.2ip is studied in [15] . The most inter- 
esting generalisation would be to models in higher dimensions. Whilst one 
cannot expect such explicit evaluations as presented here for the circle, we 
hope that the one-dimensional model clarifies some of the issues that may 
arise in higher dimensions. 
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A Appendix: Zeta-function regularisation of 
the Dirac determinant 

In this appendix, the zeta-function regularisation of det D for a self-adjoint 
operator D that is not necessarily positive definite [7] is reviewed. This 
regularisation is adapted to det(i-D), which is compared with det.D. 

A.l Positive definite D 

To begin, suppose D is a Hermitian, strictly positive operator in a finite- 
dimensional Hilbert space. The zeta-function Cd(s) of D is defined for s G C 




(4.2) 



[0,1] 1 [0,1] = 11 [0,1] 



(4.3) 
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by 

Us) = ElF> (A.1) 



where A& are the eigenvalues of D. As the eigenvalues are finitely many 
and positive, Cd(s) is well-defined and holomorphic in s. An elementary 
computation shows that 

detD = Y[\ k = e- c ^ 0) . (A.2) 



The point of flA.2j) is that the rightmost expression can be adopted as 



the definition of det D even when the Hilbert space is infinite dimensional, 
provided the spectrum of D is still discrete and the sum in (1A.1I) converges 



for sufficiently large Re s to define Cd{s) as a function that can be analytically 
continued to s — 0. The analytic continuation in s provides a prescription 
that removes the infinity from the formally divergent product Yl k A&. 

A.2 Indefinite D 

Suppose that D is an indefinite Hermitian operator in a finite-dimensional 
Hilbert space, such that the spectrum of D does not contain zero. We wish 
to express det D and det(iD) in a form similar to (IA.2|) . The new issue is to 
accommodate the negative and imaginary eigenvalues. 

Let Afc denote the eigenvalues of D, enumerated so that A^ > for k > 
and Afc < for k < 0. We define two zeta-functions by 

(dAs) = E tt + ^ S E TzxT' (A - 3a) 

fc>0 k k<0 \ k > 

q d (s) = e—i 2 E \s + eW2 E r-h^' ( A - 3b ) 

where e G {1,-1}. For integer argument these functions agree with naively 
allowing negative or imaginary eigenvalues in (lA.lj) . 
We next define the eta-function of D by 

Finally, since D 2 is positive definite, its zeta-function is defined by the re- 
placements D — > and A& — > A| in flA.ip . It follows that 

11 fc>0 fc fc<0 v /c; 
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All the functions in fl A. 3[) . (IA.4[) and (IA.5[) are well-defined and holomor- 
phic. They satisfy 



CdM = \( l + Ms/2) +1(1- ) VD (s) 
&d{s) = cos(7rs/2)C D 2(s/2) - isin(irs/2)r] D (s), 
and differentiation of (1A.6j) at s = yields 

CId(O) = ^(0) - y^(0). 



(A.6a) 
(A.6b) 

(A.7a) 
(A.7b) 



We are now ready to turn to the determinants. An elementary computa- 
tion using ( 1A.3P shows that 



E (o) 



det D = JJ A fc = e" c ^ 



det(iD) = JJ(iA fc ) = e _c « 



(o) 



(A.8a) 
(A.8b) 



Using flA.7j) . we thus have 



detD 



,"7 IB 



(0)-C D 2(0)) -|Ci a (0) 



det(iD) = e^(°) e -i^(°). 



(A.9a) 
(A.9b) 



Formulas (1A.9j) provide definitions for detD and det(i-D) that extend to 
the case when the Hilbert space is infinite-dimensional and separable and 
the spectrum of D is discrete with suitable asymptotic properties. The val- 
ues tjd{s) and Cd 2 (s/2) are in this situation defined by (IA.4j) and flA.5j) for 
sufficiently large Re s and analytically continued to s — 0. 

An important difference between det D and det(i-D) arises from the phases 
in (IA.3|) . In the definition of Qd (lA.3bj) . we chose the phases of the positive 
and negative eigenvalue terms to be the opposite for real argument, with 
the consequence that in the finite-dimensional case Qd is real- valued for real 
argument whenever the spectrum of D is invariant under D — >■ —D. In the 
definition of Cd j(E (lA.3al) . by contrast, the branch of (— l)~ s in the negative 
eigenvalue terms cannot be fixed by a similar symmetry argument, and we 
parametrised this ambiguity bye G {1,-1}, finding that e still survives in the 
final formula flA.9aj) for detD. In the finite-dimensional case rj£>(0) — Cd 2 (0) 



is an even integer, and the right-hand side of (lA.9al) is hence independent 
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of e. In the infinite-dimensional case, however, the two values of e can yield 
different regularised values for det D: an example will occur in Appendix [Bj 
The choice e = 1 is related to our regular isat ion of det (ID) since (d,i(s) = 
e ms ' 2 C l iD{s) 1 whereas the formulas in [7j make the choice e = — 1. 

B Appendix: Dirac determinant on the circle 

for U(l) 

In this appendix we evaluate det($) and det(z$) for the Dirac operator ij) 
(13.41) with the gauge group U(l), using the regularisation (IA.9I) . 

For U(l), the Dirac operator (13 .4p reduces to Ij) = i-^ + A(t), where A is a 
real- valued function of the coordinate t G [0, /], the Hilbert space is £ 2 ([0, /]), 
and both A and the domain of ]fi have periodic boundary conditions. By a 
gauge transformation we may assume A to take a constant value, which we 
denote by 2na/l, and by a further gauge transformation we may assume 
a G [0, 1). The holonomy of A is Q = e~ 2ma . Note that a is uniquely 
determined by the holonomy. 

Solving the eigenvalue equation lfiip = Xip, subject to the /-periodicity, 
shows that the eigenvalues are = 2n(k + a)/l where k G Z. We exclude the 
special case a = 0, in which one eigenvalue vanishes. We then have a G (0, 1), 
all the eigenvalues are non-vanishing, and we are in the situation covered by 
Appendix [A] 

Consider first (^2. Assuming Res > 1, we have 

= (2n/iy s (Ch(s, a) + ( H (s, 1 - a)) , (B.l) 

where the sums are absolutely convergent and the Hurwitz zeta function Ch 
is defined by [TT] 

OO j 

c*M=£(jT^r (B ' 2) 

Analytically continuing to s = and using 25.11.13, 25.11.18 and 5.5.3 
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in [TT], we find 



V(0)=0, (B.3a) 
e -^>(°) = 2s i n7ra . (B.3b) 

Consider now rjy,. Assuming again Res > 1, we have 

V#{s) = (2k/Q L | fc + ar 
k 1 1 

= (27t//)- s (Ci/(s, a) - Ch(s, 1 - a)). (B.4) 
Analytically continuing to s = and using 25.11.13 in [TT], we find 

7/0(0) = 1 -2a, (B.5) 

in agreement with [T2l §1.13]. 

Using flA.9j) with (1B.3|) and (IB.5|) . and recalling Q = e~ 2ma , we have 




det#= T" _ n ' (B-6a) 



det(i^) = 1 - Q. (B.6b) 

Note that det .0 and det(i^) depend only on the holonomy and not on 
The reason for the /-independence can be traced to the property (1B.3aj) . 



The modulus of the final result (1B.6bj) for det(i^) agrees with the cal 



culation in the physics literature of the ratio of two such determinants with 
different values of a p2]. However the phase does not agree, presumably 
due to the fact that the definition of this ratio in [13] is given as an infinite 
product that is not absolutely convergent. 

Finally, it is worth noting that when a = 0, the formulas IB.6al and IB.6bl 
all extend by continuity to give determinant equal to 0. This is consistent 
with the fact that Ip has as an eigenvalue. 
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